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Abstract 

We explore the geometry of complex networks in terms of an n-dimensional Euclidean embedding represented 
^D ' by the Moore-Penrose pseudo-inverse of the graph Laplacian (L+). The squared distance of a node i to the 

^N , origin in this n-dimensional space (/,J), yields a structural centrality index (C*{i) — 1/?,^) for node i. In turn, 

the sum of reciprocals of individual node structural centralities, J^i 1/C*(«) = J2i ^ti^ i-^- ^he trace of L+, 
yields the well-known KirchofF index (/C) , an overall structural descriptor for the network. In addition to this 
geometric interpretation, we provide alternative interpretations of the proposed indices to reveal their true 
topological characteristics: first, in terms of forced detour overheads and frequency of recurrences in random 
walks that has an interesting analogy to voltage distributions in the equivalent electrical network; and then 
as the average connectedness of i in all the bi-partitions of the graph. These interpretations respectively 
help establish the structural centrality (C*(i)) as a measure of both the overall position as well as the overall 
connectedness of node i in the network which in turn refiects the robustness of node i to random multiple 
edge failures. Through empirical evaluations using synthetic and real world networks, we demonstrate how 
the structural centrality is better able to distinguish nodes in terms of their structural roles in the network 
and, along with Kirchoff index, is appropriately sensitive to perturbations/rewirings in the network. 
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^j ' Unlike traditional studies on network robustness, that typically treat networks as combinatoric objects 

Li». . and rely primarily on classical graph-theoretic concepts (e.g. minimum cuts), we explore a geometric ap- 

proach which enables us to employ more advanced theories and techniques. More importantly, such an 
approach helps quantify and compare robustness of different nodes in the same network as well as different 
networks of comparable sizes, in terms of their local and global structures. 

1 J ■ In this work, we study a geometric embedding of networks using the Moore-Penrose (pseudo) inverse of 

rS I the graph Laplacian for the network, denoted henceforth by L+ . We show that the diagonal entries of L"*" , 

j^ ' denoted as Z^ for node i, that represent the squared distance of each node to the origin in the n — dimensional 

Euclidean space of the network embedding, provide a robust measure of a node's structural centrality, given 
as C*{i) ~ ^/iti- ^^^ particular, closer a node i is to the origin in this space, more structurally central it is. 
Moreover, the trace of L+, rr(L+) — J^i 1/C*(j), yields the Kirchoff index (/C), a structural descriptor for 
the network as a whole. Once again, lower the value of /C for a network, more compact the embedding, and 
more structurally robust the overall network is. 

In addition to the geometric interpretation of structural centrality and Kirchoff index described above, 
we provide three alternative interpretations which respectively show how C*{i) quantifies the overall position 
and connectedness of a node «, and therefore its robustness to failures in the network. We describe each of 
these separately below. 
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First, we equate structural centrality of a node i, i.e. C* {i) , to the (reciprocal of) average detour overhead 
incurred when a random walk between any source destination pair is forced to go through node i. Intuitively, 
the average overhead incurred in such detours (measured in terms of the number of steps in the random 
walk) will be lower if node i is structurally central (higher C*{i) and lower /^) and vice versa. 

Secondly, we show how C*{i) is related to electrical voltages when the network is treated as an equivalent 
electrical network (EEN). This, in turn, yields an interpretation of C*{i) in terms of the probability with 
which a random detour through i returns to the source node; referred to as the phenomenon of recurrence 
in random walk literature. More precisely, higher C*{i) implies that a random detour through node i forces 
the random walk between any source destination pair to return to the source node with lower probability, 
thereby incurring lower detour overhead. Both of these interpretations, namely average detour overhead and 
probability of recurrence, show how C*(i) quantifies the overall position of node i in the network. 

Last but not the least, we establish how C*{i) captures the overall connectedness of node i. To do so, we 
equate it to the number of nodes that i can communicate with when a sub-set of edges in the network fail in 
such a way that the network is partitioned into two connected sub-networks. A higher value of C*(z), implies 
that i is present in the larger of the two sub-networks, on an average, resulting from such partitioning. In 
other words, C*{i) reflects the immunity/vulnerability of node i towards multiple edge failures in the network, 
a distinct topological characteristic. 

Through numerical simulations using synthetic and realistic network topologies, we demonstrate that 
our new indices better characterize robustness of nodes in network, both in terms of position as well as 
connectedness, as compared to other existing metrics (e.g. node centrality measured based on degree, 
shortest paths, etc.). A rank-order of nodes in terms of their structural ccntralities helps distinguish them 
in terms of their structural roles (such as core, gateway, etc.). Also, structural centrality and Kirchoff index, 
are both appropriately sensitive to local perturbations in the network, a property not displayed by other 
centralities in literature (as shown later in this paper). 

The rest of the paper is organized as follows: we begin by providing a brief overview of several structural 
indices, characterizing node centrality as well as overall descriptors for networks, found in literature in ^ 
33] introduces a geometric embedding of the network using the eigen-space of L"*" , structural centrality and 
Kirchoff index as measures of robustness. i|4] demonstrates how structural centrality of a node reflects the 
average detour overhead in random walks through a particular node in question followed by its equivalence 
to the probability of recurrence. In ij5]we show how structural centrality captures the average connectedness 
of nodes in the bi-partitions of a network, fj^] presents comparative empirical analysis and in fJ7]the paper is 
concluded. 

2. Related Work 

Robustness of nodes to failures in complex networks is dependent on their overall position and connected- 
ness in the network. Several centralities, that characterize position and/or connectedness of nodes in complex 
networks in different ways, have therefore been proposed in literature. Perhaps the simplest of all is degree 
— the number of edges incident on a node. Degree is essentially a local measure i.e. a first oder/one-hop 
connectedness index. A second-order variant called joint- degree, given by the product of degrees of a pair of 
nodes that are connected by an edge in the network, is also in vogue. However, except in scale free networks 
that display the so called rich club connectivity |2l. llO. Ill|. neither degree nor joint-degree determine the 
overall position or the connectedness of nodes. 

A class of structural indices called betweennesses, namely shortest path/geodesic (GB) [15|,lla|, flow {FB) 
[ill and random-walk (RB) [22] betweenness respectively quantify the positions of nodes, with respect to 
source destination pairs in the network. Betweennesses reflect the role played by a node in the communication 
between other nodes and are not the measures of a node's own connectedness. 

A more sophisticated measure is geodesic closeness (GC) [la, [l6| . It is defined as the (reciprocal of) 
average shortest-path distance of a node from all other nodes in the network. Clearly, geodesic closeness is a 



p*'»_order measure of connectedness where p — {1, 2, ..., 5}, 5 being the geodesic diameter of the graph, and is 
better suited for characterizing global connectedness properties than the aforementioned indices. However, 
communication in networks is not always confined to shortest paths alone and GC being geodesic based, 
ignores other alternative paths between nodes, however competitive they might be, and thus only partially 
captures connectedness of nodes. Some all paths based counterparts of geodesic closeness include information 
centrality [26] and random- walk centrality |24| , that use random- walk based approach to measure centrality. 

Recently, subgraph centrality {SC) — the number of subgraphs of a graph that a node participates in — 
has also been proposed l9!|. In principle, a node with high subgraph centrality, should be better connected 
to other nodes in the network through redundant paths. Alas, subgraph centrality is computationally 
intractable and the proposed index in [9| approximates subgraph centrality by the sum of lengths of all 
closed walks, weighed in inverse proportions by the factorial of their lengths. This inevitably results in 
greater correlation with node degrees as each edge contributes to closed random- walks of lengths 2,4, 6, ... 
and thus introduces local connectivity bias. In a subsequent paper, Estrada et al introduce the concept of 
vibrations to measure node vulnerability in complex networks p\ . Their index called node displacement bears 
significant resemblance to information centrality and has interesting analogies to physical systems. However, 
once again the true topological significance of the centrality measure, in terms of connectedness is wanting. 

Our aim in this work, therefore, is to provide an index for robustness of nodes in complex networks, that 
effectively reflects both the position and connectedness properties of nodes and consequently, by extension, 
of the overall network. 

3. Geometric Embedding of Networks using L"*" and Structural Centrality 

In studying the geometry of networks, we first need to embed a network (e.g. represented abstractly as a 
graph) into an appropriate geometric space endowed with a metric function (mathematically, a metric space) . 
In this section we describe an n-dimensional embedding of the complex network using, the Moore-Penrose 
pseudo-inverse of the combinatorial laplacian (L+). The squared length of the position vector for a node in 
this space yields a geometric measure of centrality for the node while the sum of the squared lengths of the 
position vectors of all nodes, or the trace of L+, yields an overall robustness index for the graph. But first 
we need to introduce some basic notations. 

Given a complex network, its topology is in general represented as a (weighted) graph, G = (V,-E, W), 
where V{G) is the set of nodes representing, say, switches, routers or end systems in the network; E = {Cuv '■ 
u,v G V} is the set of edges connecting pairs of nodes representing, for example, the (physical or logical) 
communication links between the pair of nodes; and W — Wuv £ 3?+ : Cuv G E{G) is a set of weights assigned 
to each edge of the graph (here 5R+ denotes the set of nonnegative real numbers) . These weights can be used 
to represent, for example, the capacity, latency, or geographical distance, or an (administrative) routing cost 
associated with the edge (communication link) e„«. Note that if Wuv is simply or 1, we have a simple and 
unweighted graph. 

Given G = (V, E, W), we introduce ernnxn affinity matrix A = [a^] associated with G, where n — \V{G) 
is the number of nodes in G (the order of G), and atj > is some function of the weight Wij. For a simple 
graph where Wij G {0, 1}, setting a^ = Wij yields the standard adjacency matrix of the graph G. In general, 
each entry a^ captures some measure of affinity between nodes i and j: the larger Uij is, nodes i and j are 
in a sense closer or more strongly connected. Hence in general, we refer to A as an affinity matrix associated 
with G. We assume that a^ = aji, i.e. A is symmetric. For 1 < i < n, define d{i) = J^j'^ij' ^^^ refer 
to d{i) as the (generalized) degree of node i. (Note that if G is a simple unweighted graph j and A is its 
adjacency matrix, then d{i) is the degree of node i.) 

The combinatorial Laplacian of A (or the associated graph G), is defined as L = D — A, where D = 
[da] = d{i) is a diagonal matrix with d(i)'s on the diagonal. The Laplacian is a positive semidefinite matrix, 
and thus has n non- negative Eigen values A^'s. For 1 < i < n, let u^ be the corresponding eigenvector of Xi 
such that ||ui||2 — u^u^. We assume that the eigenvalues A^'s are ordered such that Ai > ... > A„ = 0. Then 
the matrix formed by the corresponding eigenvectors u^s, U = [ui, ..., u„], is orthogonal i.e. U'U = I, the 



identity matrix. More importantly, L admits an eigen decomposition L = UAU', where A is the diagonal 
matrix A = [Afi] = A^. 

Like L, its Moore-Penrose (pseudo) inverse L+ is also positive semi-definite, and admits an eigen decom- 
position of the form, L+ = UA+U', where A+ is a diagonal matrix consisting of A^^ if A^ > 0, and if 
\i = (for simplicity of notation, in the following we will use the convention A~ = if A^ = 0). Define 
X = A+^/^U'. Hence, L"'' = X'X which means that the network can be embedded into the Euclidean space 
SR" where the coordinates of node i are given by x^, the z*'* column of X. As the centroid of the position 
vectors lies at the origin in this n-dimensional space [14| . the squared distance of node i from the origin 
is exactly the corresponding diagonal entry of L"*" i.e. ||xi||2 = Z^ and the squared distance between two 
nodes i,j G V{G), \\xi — Xj||2 = Vol{G)^^Cij where Vol{G) = X]r=i '^(*) i^ called the volume of the graph 
(a constant for the graph) and Cy is called the commute time defined as the expected length of commute in 
a random walk between i and j in the network [4| • 

Based on the geometric embedding of the graph using L"*" described above, we now put forth two robust- 
ness metrics. First, a rank order for individual nodes in terms of their relative robustness called structural 
centrality, defined as: 

Definition 1. Structural centrality of node i G V{G): 

C*{^)^l/l+, yieViG) (1) 

Specifically, closer a node is to the origin in this n-dimensional space, more structurally central it is and 
consequently greater its robustness (as shown later). Next, a structural descriptor for the overall robustness 
of the network called Kirchojf index, defined as; 

Definition 2. Kirchojf index for G{V, E): 

n n 

/C(G) = Tr(L+) - 5: 4 = E VC*(z) (2) 

1=1 1=1 

Geometrically, more compact the embedding is, lower the value of IC{G) and more robust the network G is. 
It is worth pointing out that Kirchoff index has been used to model molecular strengths in the theoretical 
chemistry literature [29| , primarily because it has interesting analogies in the field of electrical networks but 
its true topological significance has scarcely been explored previously. Also, since 



n— 1 _. 

/C(G) = rr(L+) = J2- (3) 

i—l 

Kirchoff index can be thought of as a more generalized version of the much celebrated algebraic connectivity 
of the graph, measured by the second smallest eigen- value of the Laplacian i.e. A„_i [12|, |l3[. We can 
therefore use Kirchoff index to compare the robustness of two graphs with the same order and volume. 

In what follows, we demonstrate how these two metrics indeed reflect robustness of nodes and the overall 
graph respectively, flrst through rigorous mathematical analysis resulting in closed form representations and 
then with empirical evaluations over realistic network topologies. 

4. Structural Centrality, Random Walks and Electrical Voltages 

To show that structural centrality (C*(i)) indeed captures the overall position of a node, we relate it 
to the lengths of random- walks on the graph. In §4.11 we demonstrate how C*{i) is related to the average 
overhead incurred in random detours through node i as a transit vertex. Next in §4.2| we provide an electrical 
interpretation for it in terms of voltages and the probability with which a random detour through node i 
returns to the source node. 



4-.1. Detours in Random Walks 

A simple random walk {i -^ j), is a discrete stochastic process that starts at a node i, the source, visits 
other nodes in the graph G and stops on reaching the destination j [l8|. In contrast, we define a random 
detour as: 

Definition 3. Random Detour [i ^ k ^>- j): A random walk starting from a source node i, that must visit 
a transit node k, before it reaches the destination j and stops. 

Effectively, such a random detour is a combination of two simple random walks: {i — > k) followed by {k — > j). 
We quantify the difference between the random detour {i — >■ k -^ j) and the simple random walk {i — > j) in 
terms of the number of steps required to complete each of the two processes given by hitting time. 

Definition 4. Hitting Time {Hij ) ; The expected number of steps in a random walk starting at node i before 
it reaches node j for the first time. 

Clearly, Hik + Hkj is the expected number of steps in the random detour {i ^ k ^ j). Therefore, the 
overhead incurred is: 

^jjr^k^j ^ H,k + Hkj - H,, (4) 

Intuitively, more peripheral transit k is, greater the overhead in Q. The overall peripherality of k is captured 
by the following average: 

-. n n 

Alas, hitting time is not a Euclidean distance as Hij ^ Hji in general. An alternative is to use commute 
time Cij = Hij + Hji ~ Cji, a metric, instead. More importantly [20|, 

a,=VoliG){l+ + l+-l+-l+) (6) 

and in the overhead form Q, (non- metric) hitting and (metric) commute times are in fact equivalent (see 
propositions 9 — 58 in [l9j): 

^jj^^k^j = (C,fc + Cfe, - Cy)/2 = AW^''^' (7) 

We now exploit this equivalence to equate the cumulative detour overhead through transit k from ^ to Z^t^, 
in the following theorem. 

Theorem 1. 

AiJ^'^) = /+ (8) 

Proof: Using AW^^^^i = (Qk + Ckj - C,j)/2: 

-J n n 

^ ' i—l j—1 

Observing Cxy = Vol{G) (l^^ + lyy — 2/+ ) [20j and that L+ is doubly centered (all rows and columns sum 
to 0) [lj|, we obtain the proof. 

D 

Therefore, a low value of AH^'^' implies higher C*{k) and more structurally central node k is in the network. 
Theorem [T] is interesting for several reasons. First and foremost, note that: 

n 

Y, Ckj = Vol{G) {n /+ + Tr(L+)) (9) 





Figure 1; A simple graph G and its EEN. 



As rr(L+) is a constant for a given graph and an invariant with respect to the set V{G), we obtain 
^fefe '^ S?=i ^kj] lower /^^ or equivalently higher C*(A:), imphes shorter average commute times between k 
and the rest of the nodes in the graph on an average. Moreover, 



/C(G) = rr(L+) = ^;+ 



fc=l 



2nVol{G) 



fc=ij=i 



kj 



(10) 



As /C(G) reflects the average commute time between any pair of nodes in the network, it is a measure of 
overall connectedness in G. For two networks of the same order (n) and volume {Vol{G)), the one with 
lower /C(G) is better connected on an average. 

4-2. Recurrence, Voltage and Electrical Networks 

Interestingly, the detour overhead in Q is related to recurrence in random walks — the expected number 
of times a random walk (z -^ j) returns to the source i [7[. We now explore how recurrence in detours related 
to structural centrality of nodes. But first we need to introduce some terminology. 

The equivalent electrical network (EEN) 7] for G{V, E, W) is formed by replacing an edge Cy G E{G) 
with a resistance equal to w~, (see Fig. [T]) . The effective resistance [Vtij ) is defined as the voltage developed 
across a pair of terminals i and j when a unit current is injected at i and is extracted from j, or vice versa. In 
the EEN, let V^^ be the voltage of node k when a unit current is injected at i and a unit current is extracted 
from j. From [23], U]^ = d(fc)Vj!"'. Substituting k = i we get, Ul'' = d{i)Vl-' ; the expected number of times a 
random walk [i — >■ j) returns to the source i. For a finite graph G, U^'' > 0. The following theorem connects 
recurrence to the detour overhead. 



Theorem 2. 



/^Jji^k^J 



voi{G)(ut + K'-un 

d{i) 



Proof: From 27| we have, /S.H^^^^^ = d{i) ^ Vol{G) U^ . The rest of this proof follows by proving 



Ujk^mk^ljkj _jrjij^ 



From the superposition principle of electrical current, we have V^'^ — Vy" + Vy^. Therefore, 



y^k ^ ykj _ y^J 



v;'^ + (y^^ + vt'-v;f-vt) 



From the reciprocity principle, V^^ = V^^ . Therefore, V^'^ ~^Vi''~ ^l'' — ^i ■ Multiplying by d{i) on both 
sides we obtain the proof. 

D 





(a) Partition P = {S, S") (b) Spn. forest in P = {S, S') 

Figure 2; Partitions and spanning forests of a graph. 
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U, 



ki\ 



Ul'' can be interpreted as the expected extra number of times a random walk 



returns to the source i in the random detour (« ^- fc ^ j) as compared to the simple random walk {i — >■ 
j). Each instance of the random process that returns to the source, must effectively start all over again. 
Therefore, more often the walk returns to the source greater the expected number of steps required to 
complete the process and less central the transit k is, with respect to the source-destination pair (i, j). 

Therefore, Ai/^*^^ , that is the average of /\W^^^^ over all source destination pairs, tells us the average 
increase in recurrence caused by node k in random detours between any source destination pair in the 
network. Higher the increase in recurrence, i.e. AH^'^' , lower the magnitude of C*(fc) and less structurally 
central the node k is in the network. 



5. Connected Bi-Partitions of a Network 

Having established, that nodes with higher C*(i), are more structurally central in the network, we now 
turn to average connectedness of nodes in this section. To show how structural centralities of nodes capture 
their immunity/vulnerability to random failures in the network, we study their connectedness in all the 
bi-partitions of a graph. 

Definition 5. Bi-partition (P — (S, S')): A cut of the graph G which contains exactly two mutually exclusive 
and exhaustive connected subgraphs S and S' . 

Let, V{S) and V{S') be the mutually exclusive and exhaustive subsets of V{G), E{S) and E{S'), the sets 
of edges in the respective components S and S" of P and £'(5', S"), the set of edges that violate P i.e. have 
one end in 5' and the other in S". Also, let 7~(S') and T{S') be the set of spanning trees in the respective 
component sets S and S". We denote by 'P{G), the set of all bi-partitions of G{V,E). Clearly, a given 
P — {S, S') represents a state of the network in which E{S, S") have failed. A node i e V{S) stays connected 
to |V^(S')| — 1 nodes and gets disconnected from |1/^(S")| nodes. In the following relationship we show how 
structural centrality of node i is related to a weighted sum of |F(S')'| over all the bi-partitions P G P(G) of 
the network. 



Theorem 3. 



iev{S) 

z+oc Y. \ris)\\ris')\\v{s')\ 

Pev{G) 



(11) 



Proof sketch: The following result, due to Chebotarev et al. [5|, l6|, forms the basis of our proof. Let Tx 
be the set of spanning rooted forests of G{V,E) with x edges. Precisely, F^ G J-x, '^s a. spanning acyclic 
subgraph of G with the same node set as G and is composed of exactly n — x trees with one node marked as 
a root in each of these x trees. Let J-^^, be the subset of T^ in which node i is the root of the tree in which 
it belongs. Then, 

£(-^-2) - ^£(-^n-2) 



IZ = 



e(-7^«-i) 



(12) 



Here, e{-) simply represents the cardinality of the input set (see [a, Isl for details). Clearly, e{Tn~2) and 
e(J>i_i) are invariants over the set of vertices V{G) for a given graph. Hence, Z+ oc e{T^_2)- The rest of 
the proof follows from the results of the following lemmas: 

Lemma 1. Let J^^~2\p ^^ ^^^ ^^^ '^f spanning forests with n — 2 edges (or exactly two trees) rooted at node 
i in a given bi-partition P = (5*, S") and Ts, Ts' be the set of spanning trees in S and S' respectively. If 
i e V{S) then, 

eiK'-2\p) ^ \rs\ \rs'\\vis')\ 

Proof : Let Ti e Ts and T2 G Ts'- Clearly, \E{Ti)\ = \V{S)\ - 1 and \E{T2)\ = \V{S')\ - 1. As 

\ViS)\ + \ViS')\ = \V{G)\ and \E{Ti) n ^(Ta)! = 0, \E{Ti) U E{T2)\ = \V{S)\ - 1 + \V{S')\ - 1 = n - 2. 

Each such pair {Ti,T2) is a spanning forest of n — 2 edges. Given i is the root of Ti in S, we can choose 

\V{S')\ roots for T2 in 5'. There being \Ts\ \Ts'\ such pairs: el-F^^alp) = l^sl \Ts'\ \ViS')\. 

D 

By symmetry, for j G V{S') : 

£(-^^V)-|'^^ll'^^'l 1^5)1 



Lemma 2. Given V{G), the set of all bi-partitions of G: 

iev{s) 

£(-^-2)- E \r{s)\\r{s')\\vis')\ 

Pev(G) 

Proof: By definition, VF G J-n-2, F belongs to exactly one of the partitions of G. Hence, J^^*_2 = 

Upgp ^(•^n-2|p)- ^^ ^'^^ RHS is a disjoint union, counting members on both sides we obtain the proof. 

D 

Evidently, combining the results of the two lemmas above, we obtain the proof for Theorem [3l 

D 

To paraphrase, given a bi-partition P = {S,S') G V{G), such that i G V{S) and j G |F(S")|, Lemma [1] 

yields: £(-^"_2|_p)/^(-^re~2|p) ~ l^('5")l/l^('5')l- Clearly, for a given bi-partition, nodes in the larger of the 

two components of P have a lower number of spanning forests rooted at them than those in the smaller 

component and vice versa. By extension, 

it~i^ « E \ris)\\r{s'mvis')\-\vis)\) 

Pev{G) 

can be interpreted as a comparative measure of connectedness of nodes i and j. Note that for P G V{G), 
the RHS of (TTSl) is zero when nodes i and j belong to the same component of P or if |F(5')| = |F(S")| and 
positive when i G ^(5*),^ G V{S') and |V^(>5")| > |V^(5')| or vice versa. Therefore, a node i with higher 
structural centrality stays connected to a greater number of nodes on an average in a disconnected network, 
than one with lower structural centrality and is consequently more immune to random edge failures in the 
network. 

By simple extension, Kirchoff index represents the average connectedness of all the nodes when a failure of 
a subset of edges partitions the network into two halves, thereby truly reflecting overall network robustness. 
It is easy to demonstrate that of all trees of order n, the star has the lowest Kirchoff index and the root of 
the star has highest C*{i) value. Also, amongst all graphs of order n with differing volumes, the completely 
connected graph Kn has the lowest Kirchoff index. 
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(a) Abilene Topology (b) Simulated topology 

Figure 3: Abilene Network and a simulated topology. 



(c) Perturbations 



6. Empirical Evaluations 

We now empirically study the properties of structural centrality (C*(i)) and KirchofF index (we use 
K,* = K^^ henceforth to maintain higher is better). We first show in ij6.11 how a rank-order of nodes in terms 
of their structural centralities captures their structural roles in the network and then in ii6.2l demonstrate 
how it, along with KirchofF index, is appropriately sensitivity to rewirings and local perturbations in the 
network. 



6.1. Identifying Structural Roles of Nodes 

Consider the router level topology of the Abilene network (Fig. [2ta)) [1|]. At the core of this topology, is 
a ring of 11 POP's, spread across mainland US, through which several networks interconnect. Clearly, the 
connectedness of such a network is dependent heavily on the low degree nodes on the ring. For illustration, 
we mimic the Abilene topology, with a simulated network (Fig. [Sfb)) which has a 4-node core {ui, ...,U4} 
that connects 10 networks through gateway nodes {ws, ...,wi4} (Fig|3ljb)). 

Fig. |4] shows the (max- normalized) values of geodesic closeness (GC), subgraph centrality (SC) and 
structural centrality C* for the core {ui, ..., U4}, gateway {us, ..., U14} and other nodes {wis, ..., ues} in topology 
(FigElJb)). Notice that v^ and ug, two of the gateway nodes in the topology, have the highest values of degree 
in the network i.e. {d{v^) = d{vQ) = 10) while W14 has the highest subgraph centrality (SC). In contrast, C*{i) 
ranks the four core nodes higher than all the gateway nodes with vi at the top. The relative peripherality 
of W5 , vq and vn as compared to the core nodes requires no elaboration. As far as geodesic centrality (GC) 
is concerned, it ranks all the nodes in the subnetwork abstracted by W5, namely W15 — W23, as equals even 
though U22 and U23 have redundant connectivity to the network through each other and are, ever so slightly, 
better connected than the others - a property reflected in their C*{i) rankings. 

We see similar characterization of structural roles of nodes in two real world networks in terms of structural 
centrality: the western states power-grid network [28| and a social network of co-authorships [23|, as shown 
through a color scheme based on C*{i) values in Fig. [S] Core- nodes connecting different sub-communities 
of nodes in both these real world networks are recognized effectively by structural centrality as being more 
central (Red end of the spectrum) than several higher degree peripheral nodes. 
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Figure 4: Max-normalized centralities for simulated topology. 



6.2. Sensitivity to Local Perturbations 

An important property of centrality measures is their sensitivity to perturbations in network structure. 
Traditionally, structural properties in real world networks have been equated to average statistical properties 
like power-law/scale- free degree distributions and rich club connectivity [3, [lOl, lUl • However, the same degree 
sequence D = {d{l) > d{2) > ... > d{n)}, can result in graphs of significantly varying topologies. Let G{D) 
be the set of all connected graphs with scaling sequence D. The generalized Randic index Ri{G) |3ll25|: 



RiiG)= J2 ^W^(J') 



(13) 



eijeE(G) 



where G G G{D), is considered to be a measure of overall connectedness of G. Higher Ri{G) suggests 
that nodes of higher degrees connect with each other with high probability thereby displaying the so- 
called rich club connectivity (RCC) in G |21| . Similarly, the average of each centrality/betweenness index 
{GC,SC,GB,RB averaged over the set of nodes), is in itself a global structural descriptor for the graph 
G [9| . We now examine the sensitivity of each index with respect to local perturbations in the subnetwork 
abstracted by the core node vi and its two gateway neighbors v^ and vq. 

First, we rewire edges 615,5 and ee.i to 6154 and 66,5 respectively (PERT-I Fig. Eljc)). PERT-I is a degree 
preserving rewiring which only alters local connectivities i.e. neither individual node degrees nor average node 
degree changes. Fig. M,s.) and Fig. [7] respectively show the altered values of centralities {C*,GC,SG) and 
betweennesses, geodesic and random- walk i.e. {GB,RB), after PERT-I. Note, after PERT-I, W15 is directly 
connected to vi which makes C*{vi^) comparable to other gateway nodes while SC{vi5),GB{vi^),RB(vi^) 
seem to be entirely unaffected. Moreover, PERT-I also results in wg losing its direct link to the core, reflected 
in the decrease in C*{vq) and a corresponding increase in C*{vq). C*{i), however, still ranks the core nodes 
higher than U5 (whereas SC, GB, RB do not) because PERT-I being a local perturbation should not affect 
nodes outside the sub-network — vi continues to abstract the same sub-networks from the rest of the 
topology. We, therefore, observe that C*{i) is appropriately sensitive to the changes in connectedness of 
nodes in the event of local perturbations. But what about the network on a whole? 

Let G and Gi be the topologies before and after PERT-I. Gi is less well connected overall than G as 
the failure of 65^ in Gi disconnects 19 nodes from the rest of the network as compared to 10 nodes in G. 
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(a) The western-states power grid network [28] (b) A network of co-authorships in network sciences [231 

Figure 5: Real world networks: Red — > Turqoise in order of decreasing C*(i). 
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Figure 6: Max-normalized values of structural centrality, geodesic closeness and subgraph centrality for core, gateway and some 
other nodes. 



However, 



Ai?i(G ^ Gi) = ^i(gi^) 73(g) = 0.029 



i?i(G) 

as the two highest degree nodes {v^ and vq) are directly connected in Gi (see TABLE [1] for the sensitivity 
of other centrahty based global structural descriptors). In contrast, AIC*{G -^ Gi) ~ —0.045, which rightly 
reflects the depreciation in overall connectedness after PERT-I (recall /C*(G) — K.^^{G). Table [1] shows the 
changes in the average of all centrality and betweenness indices post PERT-I. 

A subsequent degree preserving perturbation PERT-II of Gi, rewiring 622.23 and 624,25 to 622.25 and 
623,24, to obtain G2, creates two cycles in G2 that safeguard against the failure of edge 65.6- This significantly 
improves local connectivities in the sub-network. However, Ai?i(Gi — >■ G2) = (and average SC decreases) 
while A/C*(Gi — s- G2) = 0.036 which once again shows the efficacy of Kirchoff index as a measure of global 
connectedness of networks. 
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Figure 7: After PERT-I: Max-normalized values of structural centrality, geodesic and random-walk betweennesses for core, 
gateway and some other nodes. 

Table 1: Sensitivity to local perturbations, X = l/n^"X(j): Avg. node centrality for a network. 



# 


Structural Descriptor 


PERT-I 


PERT-II 
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K*{G) (orC*) 
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Ri{G) 
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3 


GC 
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4 


SG, GB 
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5 


RB 


t 


t 



7. Conclusion and Future Work 

In this work, we presented a geometric perspective on robustness in complex networks in terms of the 
Moore-Penrose pseudo-inverse of the graph Laplacian. We proposed structural centrality (C*(i)) and Kirchoff 
index {IC{G)) that respectively reflect the length of the position vector for a node and the overall volume 
of the graph embedding and therefore are suitable geometric measures of robustness of individual nodes 
and the overall network. Additionally, we provided interpretations for these indices in terms of the overhead 
incurred in random detours through a node in question as well as in terms of the recurrence probabilities and 
voltage distribution in the EEN corresponding to the network. Both indices reflect the global connectedness 
properties of individual nodes and the network on a whole, particularly in the event of multiple edge failures 
that may render the network disconnected. Through numerical analysis on simulated and real world networks, 
we demonstrated that C*{i) captures structural roles played by nodes in networks and, along with Kirchoff 
index, is suitably sensitive to perturbations/rewirings in the network. In terms of computational complexity, 
structural centrality compares well with other geodesic and all-paths based indices in literature (see Table [5]) 
and performs better than random- walk betweenness in the asymptotic case. In future, we aim at investigating 
similar metrics for the case of strongly connected weighted directed graphs to further generalize our work. 
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Table 2; Taxonomy and computational complexities of centrality measures (for all nodes). 



# 


Measure 


Paths covered 


Complexity 


1. 


Degree 


- 


0{m) 


2. 


GC,GB 


Geodesic paths 


0{n') 


3. 


C* 


All paths 


0{n') 


4. 


sc 


AU paths 


O(n^) 


5. 


RWB 


AU paths 


0{m + n)n^ 


6. 


FB 


AU paths 


O(m-'n) 
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